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Abstract 



We analyze a class of supersymmetric SO(10) grand unified theories with type II seesaw 



for neutrino masses, where the contribution to PMNS matrix from the neutrino sector has an 
exact tri-bi-maximal (TBM) form, dictated by a broken S4 symmetry. The Higgs fields that 



determine the fermion masses are two 10 fields and one 126 field, with the latter simultaneously 
contributing to neutrino as well as charged fermion masses. Fitting charged fermion masses and 
the CKM mixings lead to corrections to the TBM mixing that determine the final PMNS matrix 
with the predictions #13 ~ 4° — 6° and the Dirac CP phase to be between —10° and +15°. We 
also show correlations between various mixing angles which can be used to test the model. 



I. INTRODUCTION 

Understanding neutrino masses and mixings is an integral part of our attempts 
to unravel the flavor puzzle in particle physics. During the past decade, the large amount 
of information on neutrino masses and mixings gained from the study of accelerator, 
reactor, solar, and cosmic ray neutrino observations has given a strong forward momentum 
to this journey. Several crucial pieces of the puzzle must still be found before we can begin 
to have a complete picture at hand; among them are the nature of the neutrino masses 
(Dirac vs Majorana), the mass hierarchy (normal vs inverted), the mixing angle #13, and 
the CP phases. 

Of the large number of new experiments that are under way to answer these 
questions, the T2K experiment has recently announced a possible indication of a non- 
zero value for # 13 l| , which has caused a great deal of excitement in the field. The T2K 
lower limit, if correct, is not far below the current experimental bound from the CHOOZ 
experiment [2j and has important theoretical implications. The MINOS experiment has 
also seen an excess of electron events which could be indicating a non-zero #13 [3i], and 
their allowed range for #13 overlaps with the T2K one. There have also been analyses 

j. Additionally, other experiments 



5l, anc. 
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several recent papers have 



of existing oscillation data suggesting a non-zero #13 [ 

are currently searching for this important parameter 

attempted to explain the T2K values within different models [6j] ; there is hope the situation 

will become much clearer in near future. 

A non-zero value for #13 has profound implications for our understanding of 
the physics of neutrino mass. It is, for example, well known that maximal atmospheric 
neutrino mixing (tan 6*23 = 1) suggests an underlying discrete \i — r symmetry (denoted 
by Z 2 ^- T ) in the neutrino mass matrix, which, when exact, leads to vanishing 9 13 [3] . 
Depending how this symmetry is broken (e.g. in the //-sector or e-sector), the resulting 
value of 9 is can either be very small or not so small. The neutrino mass matrix has 
also been suspected to have a larger symmetry beyond this from the observation that 
the current values of the solar mixing angle seems to have a geometric value (tan 12 = 
l/v2). The resulting lepton mixing (PMNS) matrix is known as the tri-bi-maximal 
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mixing matrix [8] (TBM for short). This symmetry is often denoted by two Z 2 symmetries 



Z2.S x ^2,/i-r [9]. This full symmetry leads to zero #13 and restricts the form of the neutrino 
mass matrix (to be called TBM matrix) to 



/, 



M v 



\ 



V 



(i) 



be c 
c a + b c — a 
c c — a a + b I 

which is given by only three parameters. In fact, in the above matrix, one could set b = 0, 
without changing the TBM PMNS matrix. It only affects the masses of the neutrinos. 
This matrix is very different from the known mass matrices in the quark sector and could 
be a possible clue to a unified understanding of the quark-lepton flavor puzzle. A non- 
zero #13 suggests that the TBM PMNS mixing is not precisely the right form, and that 
"large" corrections to both the /i — r symmetry and TBM matrix must be present; these 
factors could eventually guide us towards a complete determination of the neutrino mass 
matrix. Once this is accomplished, we will have passed a major milestone in uncovering 
the physics of neutrino mass and possible underlying symmetries of the lepton sector. Of 
course, if observations require that the corrections to TBM mass matrix are "large", it 
would not be too implausible that the symmetries described above may only be illusory 
and some other mechanisms may be at work. 

To explore what other scenarios could lead us to the desired neutrino mixings 
with a "large" #13, recall that there is a large class of predictive SO (1 0) g rand unified 



models 
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j-ll3J] in which neutrino masses arise out of a type II seesaw 14j mechanism. 
These models provide a natural way to understand a large atmospheric mixing angle 
not from some symmetry, but rather from the dynamical property that in grand unified 
theories, the bottom and tau masses become nearly equal at GUT scale [11] . When these 
models are analyzed for the full three generation case, one finds, in addition to a large 
#12, that #13 is also generally "large" [12] . Though the first of these results were obtained 
without quark CP violation, these models have since been studied in much greater detail 
and including the phase. These full CP-violating models do confirm the above results 



including a "large" #13 as well 15], but at the cost of severely restricting the parameter 
space. It turns out, however, that a slight extension of the Higgs sector by the addition of a 
120 Higgs multiplet 16] considerably broadens the parameter space while still preserving 



the "large" 6*13 prediction. 

Recently, an interesting connection to the standard TBM model discussion 
has been noted: the type II seesaw formula for neutrino masses allows a TBM form 
for the neutrino mass matrix by simply a choice of fermion basis, with no additional 
symmetries 17] : corrections to the TBM form then arise from the form of the charged 
lepton matrix, which, in our case, is determined by the 50(10) constraints from quark 
masses and mixingso Strictly speaking, no bottom-tau unification is invoked in this 
approach. Detailed numerical analyses of these models have been carried out and lead 
to excellent fits for models with 10, 126 and 120 Higgs fields 17], ll9|, and, yet again, 
a large 6*13 is predicted. One could therefore construe the "large" #13 prediction of these 
models as an indication of grand unified origin of neutrino masses (especially of the kind 
noted), which was anyway suspected as a possibility due to the near-GUT seesaw scale. 

As mentioned above, the near-tri-bi-maximal PMNS form in this class of GUT 
theories is related to the dynamics of the model rather than to any symmetry. Of course, 
to understand the particular Yukawa textures, one may need to invoke some symmetries, 
but still those symmetries are not directly related to the #13 value. We are therefore 
faced with two contrasting but attractive approaches to current neutrino observations: 
one based on leptonic symmetries, and another based on grand unification hypothesis. 
It is clearly important that more work be done to uncover which is the path chosen by 
nature; in this paper, we further investigate the grand unification approach. 

One straightforward way to establish that a "large" #13 is a generic prediction 
of 50(10) models with type II seesaw and their associated dynamical properties, rather 
than a symmetry, is to study more of such models and establish their predictions. A 
particularly simple class of models are defined by the minimal choice of Higgs fields 10 
and 126, together with either a 120 16[ (as already noted) or an extra 10 J20j contributing 
to fermion masses. The latter class of models, to the best of our knowledge, has not been 
thoroughly scrutinized numerically. In this paper, we focus on them, since, as has been 
recently pointed out [20j], they seem to give qualitatively the right picture for not just 
neutrino masses but quarks as well. It was shown in Ref. 20j that reasonably well-satisfied 
1 For charged lepton corrections to tri-bi-maximal mixing outside the framework of GUT theories, see 
Ref. [18J. 



versions of the GUT scale relations m& ~ m T and m M ~ — 3m s emerge out of an S4 flavor 
symmetry in an SO (10) GUT model of the above type. It was also noted in this paper that 
the TBM form of the neutrino mass matrix is dictated by the S4 symmetry breaking. From 
our quantitative analysis of this model, we first find that the Yukawa texture predicted 



by the minimal version of the model [20] needs to be supplemented by additional effective 
GUT scale Yukawa couplings in order to come close to observations. The improved model 
has only twelve parameters and is therefore predictive in the neutrino sector. We find 
that the model leads to a prediction for (9 13 ~ 4° — 6° and Dirac CP phase between —10° 
and +15°; this value of #13 supports the generic expectation for this class of theories, as 
was anticipated above. We also argue that there is a definite kind of correlation between 
the #13 and #23; which can be different from non-GUT symmetry-based approaches to 
#13. It is also interesting that this value is consistent with the recent T2K range for this 
parameter. 

II. DETAILS OF THE MODEL 

The class of 5*0(10) models in which we are interested have two 10 Higgs fields 



(denoted by H, H') and a pair of 126 + 126 (denoted by A + A). The SO (10) invariant 
Yukawa couplings of the model are given by: 

C Y = h ipipH + titfpipH' + /oV^A (2) 

where ?/>'s denote the 16 dimensional spinors of SO (10) which contain all the matter 
fields of each generation, and so there are three such fields, though we have suppressed 
the generation indices. The Yukawa couplings are 3x3 matrices in generation space. The 
effective Yukawa couplings fo,ho, h' are assumed to have descended from a higher scale 
theory which has an S4 symmetry broken by flavon fields with particularly aligned vacuum 



expectation values (vevs) (see e.g. Ref. (20J|). We do not need to know the detailed form 
of these flavon interactions for our analysis in this paper, and we will simply write down 
the effective form of the h, h ; , f that follow from it. Before doing that, we wish to point 
out that it is the /o coupling which is responsible for neutrino masses via type II seesaw 
mechanism, and it also contributes to charged fermion masses. We can choose it to give 



the tri-bi-maximal form for M. v either by a choice of the basis of matter fields [17], |20] 



/o ocM 



/-,; 



(3) 



/ 



or by the breaking of the S4 symmetry [20| or other symmetry group 21] . This puts 
the neutrino mass matrix in a form that, upon diagonalization, leads to tri-bi-maximal 
mixing prior to charged lepton corrections. This requires that we have the /o coupling in 
the form: 

' m 1 mi ' 

nil m o t^i ~ m o 
\ nil m \ — m o m o 

where we have used Eq. (pQ) with rescaled variables a = kttiq and c = mxi\. Note that 
as in Ref. [20|, we have taken b = in Eq. ([ID. The proportionality constant between 
/o and ./VC is determined by the left triplet vev in 126 responsible for type II seesaw. 
Note that Upmns = V^V u (where Ve and V v are the unitary matrices that diagonalize the 
charged lepton and neutrino mass matrices respectively) so that we will necessarily have 
corrections to the TBM mixing coming from the charged lepton mass matrix. Note further 
that since the /o matrix also contributes to the quark and charged lepton masses, neutrino 
masses and quark masses are connected, making the model predictive. The formulae for 
the quark and charged lepton masses in this model are given by: 



(4) 



M u ~- 


= h + r 2 f + r 3 ti, 


M d ~- 


tanp 


M t = 


= -H-(h-3f + h') 



where /, h, h! are related to / , ho, h' through Higgs vevs [20|. In Ref. [20(], the ^4 symmetry 
constrains h to be a rank one matrix of the form: 

'0 N 
h = , (5) 

V My 



and the form of h' to be 



u 



1 8 -(J^ 

5 
v -5 y 



(6) 
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The parameters mo, mi in Eq. ([3]) are chosen to be real. The parameters ri t2 ,3 in Eqs. (jlj) 
represent the ratio of the different standard model (SM) doublet vevs in the theory. There 
are three SM doublets and hence six vevs; three of them are absorbed to redefine the 
Yukawa couplings from dimensionless ho,h' , fo to h,hf, f with dimensions of mass. Since 
we have chosen h to be in the form above due to the S4 symmetry breaking, it has only 
one parameter. The / matrix has two real parameters, and h! has only one parameter, 
which can be chosen to be complex, for a total of eight parameters in the char ged fermion 
sector. While this model has a number of attractive features as noted in Ref . 20| , it fails 
to reproduce some details of the quark mixings, e.g. both V c b and V u b come out to be 
much too small compared to their extrapolated values at the GUT scale for all tan/3; the 
CKM phase also comes out too small. We therefore amend this model by adding extra 
structure to the h! matrix while keeping all other couplings as they were. We choose h! 
to have the form: 

I 5" 5 -5 + 5'\ 
5 d 

-5 + 5' d 



ti 



\ 



(7) 



/ 



Mo 



M, 



n 

tan/3 

tan/3 
/ 



which can be generated by choice of flavon fields and the alignment of their vevs. The 
neutrino mass matrix is unaffected by this addition, but the quark and charged lepton 
mass matrices are now 

( 5" -3m! + 5 -3m! - 5 + 5' N 

— 3m! + 5 — 3m 3m — 3mi + d 

t — 3m-i — 5 + 5' 3m — 3mi + d — 3m + M j 

( 5" mi + 5 m-i-5 + 5' \ 

mi + 5 mo —mo + mi + d 

\ mi — 5 + 5' —mo + mi + d mo + M < 

r^5" vimi + r%5 r^viii — rz(5 — 5') 

r 2 mi + r 3 5 r 2 m -r 2 m + r 3 (mi + d) 

\ r 2 mi - r 3 (5 - 5') -r 2 m + r 3 (mi + d) r 2 m + M j 



(8) 



M u 



III. PREDICTIONS OF THE MODEL 



The model has eleven parameters if we choose all except 5' real (twelve pa- 
rameters when we allow 5 complex to study the allowed range of Dirac CP phase). Recall 



13| . In that 



that the model with 10, 126 and 120 has a total of seventeen parameters 
sense ours is a more economical one and is quite predictive. Before proceeding with the 
numerical analysis discussion, we note a few results that can be derived analytically if we 
assume the hierarchy M 3> tuq, d 3> mi, 8, 8', 8": 



i'i 



m a 



tan/3 
r 2 m , r 2 



^; M 

m t 

m b m c 



m t : 



m 



m t m s 



(9) 



m t m s irif, 

Diagonalizing the matrices in Eqs. (JHJ) gives the charged fermion masses, and the combi- 
nation V^Vd (where V u and V& diagonalize the up- and down-sector quark masses respec- 
tively) gives the CKM matrix. Approximate expressions for the mass eigenvalues and the 
CKM mixing matrices are given by 



M, 



D 



n 



tan/3 



ai _|_ (-3toi+5) 2 
"t" <5"+3m 



— 3m r 



(Srrn+Sf 
3m +5" 



[3(mo-mi)-i 
M 



M°~ 



tan/3 



5" + 



(mi+i) 2 



-3m + M j 
\ 



-mo 



'«o + £# 



(-toq+toi+cQ 2 
M 



(10) 



m + M J 



M, 



D 



r 3 8" + (r2 A 7 +r ^ )2 

° r^o" — r-2,mo 



■ (r2mi+r3<5) 2 

roTUn + v ,, 

* u r2mo— r^b" 



-r2ino+r3(mi+d)]' 2 
M 



r 2 m + M j 



and 



V d 



1 

m\+& 
5" -mo 

mx — S+S' 
\ 5"-m -M 

( 1 



mi+S mi— 5+8' i 



mo- 5" 

1 

mp— mi- rf 
M 



-5"+m +M 

-mp+mi+rf 
M 



K 



r2mo—r38" 



rgmi+rgjj -i 

rsS"— r^rno 

r2»»i-r3(i-i') rgmp— T3(mi+d) 
\ r 3 8"-r 2 m -M M 



J 

r2mi—r-z(8—5') \ 
— r3<5"+r2>no+M 
-r2m +r3(mi+d) 
M 



1 



:n^ 



/ 



Additionally, note that the resulting corresponding expression for the Cabibbo angle is 

m 1 + 5 r 2 mi + r 3 5 



V,, 



(12) 



rriQ — 5" r2mo — r^5" 
Using a sufficient set of the individual expressions for Vy and rrif above, as 
well as the ratio mi/ mo from the neutrino sector (to be discussed later), we solve a 
system of equations against experimental values for the charged fermion masses and quark 
mixings to find an analytical solution with approximate values for the input parameters; 
this solution is then used to generate predictions for the neutrino sector and is made 
statistically robust through numerical analysis. A best fit value for the input parameters 
is given in Table [J and the resulting mass and mixing parameter values are given in 
Table [TT1 Note that while we get a higher value for ra& and slightly lower values for m s 





tan/3 = 10 


tan /3 = 55 


M (GeV) 


88.2 


106.2 


m (GeV) 


1.435 


1.382 


rrn (GeV) 


0.275 


0.275 


5 (GeV) 


0.2850 


0.2605 


5' (GeV) 


0.463 - 0.279* 


0.529 - 0.335z 


5" (GeV) 


-0.0652 


-0.0767 


d (GeV) 


3.78 


4.31 


7*1 / tan (3 


0.0153 


0.0159 


T2 


0.130 


0.129 


n 


-0.06 


-0.07 



TABLE I. Best fit values for the model parameters. Note that adding a small imaginary part 
to S will give us a non-negligible Dirac CP phase, as shown later in Figured) 



and m c , a 
literature 



1 are within reasonable statistical deviation from extrapolated values in the 



22]. Similarly, our predictions for ma and m u are somewhat higher than those 



obtained in Ref. 22], but we believe there could easily be instanton corrections to the light 
quark masses, which could change these extrapolated values. It is nevertheless remarkable 
that we are able to reproduce all other parameters in the charged fermion sector so well. 





tan /3 = 10 


tan (3 = 55 


best fit 


RG extrapolated 


best fit 


RG extrapolated 


m e (MeV) 


0.3587 


n Qt;oc+0.0003 


0.3563 


S565+ - 0002 


m^ (MeV) 


75.6865 


75 671 ^+0-0578 


75.3359 


75.2938+.°,;°^ 


m T (GeV) 


1.2927 


1 9QOO+0.0013 

i.zyzz_ 00012 


1.6272 


1 R9Q9+0- 0443 

i.ozyz_o.02 94 


m d (MeV) 


3.8587 


1 5036+ - 4235 

J..OUOD_ 2304 


4.0202 


1.4967+^ 


m s (MeV) 


23.6026 


29.9454+t300i 


23.1619 


29.oloO_4 4gg7 


m b (GeV) 


1.3726 


1.0636_£S£_ 


1.7078 


1 4167+ - 4803 


m u (MeV) 


1.9772 


0.7238+|];^5 


3.6311 


0.7244+S 


m c (MeV) 


177.3862 


210.3273+2?;°|! 


177.6719 


210.5048ig : }gg 


m t (GeV) 


88.3886 


on /(qqq+30.2676 

a/.Asaa_ li 7686 


106.3806 


nc: i 4 o fi +69.2836 
95.148b_2Q.659o 


V_ fl 


0.2230 


0.2243 ± 0.0016 


0.2233 


0.2243 ± 0.0016 


K6 


0.0032 


0.0032 ± 0.0005 


0.0032 


0.0032 ± 0.0005 


V cb 


0.0349 


0.0351 ± 0.0013 


0.0352 


0.0351 ± 0.0013 


<5CKM 


-64.35° 


(-60 ±14)° 


-61.84° 


(-60 ±14)° 



TABLE II. The best fit values of the quark and charged lepton masses and the most relevant 
quark mixing parameters. The ltr experimental values extrapolated by MSSM renormalization 
group (RG) equations to the GUT scale [22j are also shown for comparison. 



For the neutrino sector, the structure of the mass matrix in the model of 
Ref. [_0| and in this amended model is given by 

/ 



M„ 



K 



mi mi 

nil mo mi — nio 

,2 



\ 



(13) 



\ nil mi — mo mo cos a I 
where ft is a scaling factor determined from experimental data, and a is the mixing angle 
for the third generation matter fermion ip with the vector-like field ipy specific to the 
model. The limit a = gives the strict TBM form for the neutrino mass matrix given in 
Eq. ([_} when the mass eigenvalues are 



m,. 



— ftmi, 



m 



I '2 



2ftmi, 



m 



"3 



ft (2mo — mi) 



(14) 
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and the solar-to-atmospheric mass-squared ratio is given by 






3 A 2 



41-A 



(15) 



where A = rrii/mo- To fit the experimental data, Amg/Am^ ~ 0.03, which corresponds 
to A ~ 0.2 from Eq. (fTBTl . This constraint can be relaxed for the case of a 4 1 0, which is 
already required for the top Yukawa coupling (ex sin 2 a) to be non-zero |20|, though we 
do not have much freedom for the value of A anyway, as it is tightly constrained by the 
quark sector. Numerically, we find that the allowed range of a is 5° — 25° in order to fit 
the observed neutrino data. 

Noting the charged lepton rotation matrix from the ansatz given by Eq. (jSj): 

/ i 3mi- <5 

5"+3m -5"-3mo+M 

(16) 



V, 



1 

-3'mi+<5 



5"+3m 
1 



<5"+3m 
, —3mi — S+S' 3(mi—mo)—d -< 

\ 5" " " 



Smx-S+S' \ 
-S"-3m +M 

3(mo—mi)+d 
M 



/ 



'+3m -M M 

and given the TBM form of the matrix that diagonalizing the neutrino mass matrix: 



Vn 



TUB 



( 



\ 







\ 



/ 



(17) 



we can write an approximate analytical form of the PMNS neutrino mixing matrix: 



Upmns — Ve Vtbm — Vtbm + 

/ J_ ( 3m 1 -5 _ 3^ml+^-^ /, ' 

y/6 \5"+3m 

1 f 3(m 1 -m )-d _ 2(3toi+<5-<5 / *) 
V VE \ M M-3m -8" 



M-3m -8" 
'3(mo- m\)+d , 2(3rrti — S) 
M ' 8"+3m 



3m 1 ±S-5' 



1 I -3mi 

v^ \S"+3m "T" M-3m -6" 

1 /3(mi- mo) — d 

73 V 



M 
1 f 3(m —mi)+d 



I 3rrti— <5 \ 
"^ 8"+3rn ) 

3mi+(5-<5 / * 
M-3mo-5" 



1 / 3mi-S _i_ 3m 1 +5-5 r 
^/2 V<5"+3m ^ M-3m -6" 
3(mi—m())—d 

73m 

3(mi—mo)—d 
s/2M 



The exact numerical results for neutrino mixing corresponding to the quark 
sector fit from above are given in Figures [T]|U Figures [1] and [2] show the relationships 
between (6> 2 3,6> 13 ) and (9 12 , 13 ), respectively. Note that the value of #i 3 is large, though 
not so large as the 6°-8° central value of T2K result gives. Also note that the atmospheric 
mixing angle 823 is always larger than the maximal value of 45°; this agrees with the 
analytical form given in Eq. (fi~8l) in which the corrections from the charged lepton sector 



11 




FIG. 1. Correlation between #13 and #23 predicted in our model satisfying all the charged fermion 
sector constraints. The thin (thick) dotted vertical line is the current 2a (3<r) upper limit for 
the atmospheric mixing angle. 
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34 35 
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FIG. 2. Correlation between #13 and 612 predicted in our model. The solid and dotted vertical 
lines are the current la and 3a limits respectively for the solar mixing angle. 



are always positive. Figure [3] shows the correlation between (6*13, Amg/Am^j. The solid 
and dotted lines are the current la and 3a limits for the best fit values of the observed 



neutrino oscillation parameters 23( (using the new reactor neutrino fluxes). Figure H] 
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FIG. 3. Correlation between Am^/Am^ and #13 predicted in our model. The solid and dotted 
horizontal lines are the current la and 3<r limits respectively for the solar-to-atmospheric mass 
squared ratio. 
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FIG. 4. Correlation between the Dirac and CKM CP phase. The dotted vertical line is the 
central value of the observed CKM CP phase. The spread in phase in the model arises from 
complexifying the parameter S so that it does not pull the CKM phase out of its allowed range. 

shows the correlation between the Dirac and CKM CP phases. These correlations between 
different mixing parameters could be used to test the model once the current uncertainties 
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in both 9 12 and #23 are reduced and a more precise value for #13 has been determined. 
Our model also predicts small Majorana phases (~ 1°). 

IV. SUMMARY 

We have analyzed the predictions of an 50(10) model with type II seesaw for 
neutrino masses and Yukawa couplings involving two 10 Higgs fields and one 126 Higgs 
field, with all the couplings derivable from a broken S4, symmetry. The model has at most 
twelve parameters and is thus a relatively economical one when compared to other models 
discussed in the literature. It gives a fairly good fit to the charged fermion masses as well 
as an excellent fit to the CKM parameters, and it also predicts the neutrino mixing angles 
$12, #23 as well as Am^/Am^ in agreement with observation. Furthermore, it predicts 
a value for # 13 between 4° — 6°, near the lower end of the current T2K allowed range. 
With more accurate determination of #13 - especially its correlation with #23, which our 
model predicts to be strictly larger than 45° - the model could be tested. Finally, the 
model predicts a normal hierarchy for the neutrinos and hence an effective neutrino mass 
in neutrino-less double beta decay, which is a few milli-electron-volts and is thus not 
observable in the current round of the searches for this process. 
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